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■ We present in this work the study of the linear perturbations of the 2 + 1-dimensional circularly 

symmetric solution, obtained in a previous work, with kinematic self-similarity of the second kind. 
We have obtained an exact solution for the perturbation equations and the possible perturbation 
C^l ' modes. We have shown that the background solution is a stable solution. 

>v 

I. INTRODUCTION 

One of the most outstanding problems in gravitation theory is the study of the relation that exists between the 
critical phenomena and the process of black hole formation. The studies of non-linearity of the Einstein field equations 
near the threshold of black hole formation reveal very rich phenomena [1] , which are quite similar to critical phenomena 
in Statistical Mechanics and Quantum Field Theory [2]. In particular, by numerically studying the gravitational 
' collapse of a massless scalar field in 3 + 1-dimensional spherically symmetric spacetimes, Choptuik found that the 
■"sj" , mass of such formed black holes takes a scaling form, 
O ■ 

>rv M BH = c{p){p-p*)\ (i) 

o : 

' where C (p) is a constant and depends on the initial data, and p parameterises a family of initial data in such a way 
that when p > p* black holes are formed, and when p < p* no black holes are formed. It was shown that, in contrast to 
C(p), the exponent 7 is universal to all the families of initial data studied. Numerically it was determined as 7 ~ 0.37. 
- J The solution with p = p* , usually called the critical solution, is found also universal. Moreover, for the massless 
!— 1 , scalar field it is periodic, too. Universality of the critical solution and exponent, as well as the power-law scaling 
of the black hole mass all have given rise to the name Critical Phenomena in Gravitational Collapse. Choptuik's 
J> \ studies were soon generalised to other matter fields [3,4], and now the following seems clear: (a) There are two types 
of critical collapse, depending on whether the black hole mass takes the scaling form (1) or not. When it takes the 
scaling form, the corresponding collapse is called Type II collapse, and when it does not it is called Type I collapse. 
In the type II collapse, all the critical solutions found so far have either discrete self-similarity (DSS) or homothetic 
self-similarity (HSS), depending on the matter fields. In the type I collapse, the critical solutions have neither DSS nor 
HSS. For certain matter fields, these two types of collapse can co-exist, (b) For Type II collapse, the corresponding 
exponent is universal only with respect to certain matter fields. Usually, different matter fields have different critical 
solutions and, in the sequel, different exponents. But for a given matter field the critical solution and the exponent 
are universal. So far, the studies have been mainly restricted to spherically symmetric case and their non-spherical 
linear perturbations. Therefore, it is not really clear whether or not the critical solution and exponent are universal 
with respect to different symmetries of the spacetimes [5,6]. (c) A critical solution for both of the two types has one 
and only one unstable mode. This now is considered as one of the main criteria for a solution to be critical, (d) The 
universality of the exponent is closely related to the last property. In fact, using dimensional analysis [7] one can show 
that 

where k denotes the unstable mode. 
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From the above, one can see that to study (Type II) critical collapse, one may first find some particular solutions 
by imposing certain symmetries, such as, DSS or HSS. Usually this considerably simplifies the problem. For example, 
in the spherically symmetric case, by imposing HSS symmetry the Einstein field equations can be reduced from 
PDE's to ODE's. Once the particular solutions are known, one can study their linear perturbations and find out 
the spectrum of the corresponding cigenmodes. If a solution has one and only one unstable mode, by definition we 
may consider it as a critical solution (See also the discussions given in [8] ) . The studies of critical collapse have been 
mainly numerical so far, and analytical ones are still highly hindered by the complexity of the problem, even after 
imposing some symmetries. 

Lately, Pretorius and Choptuik (PC) [9] studied gravitational collapse of a massless scalar field in an anti-de Sitter 
background in 2 + 1-dimcnsional spacetimes with circular symmetry, and found that the collapse exhibits critical 
phenomena and the mass of such formed black holes takes the scaling form of Eq.(l) with 7 = 1.2 ± 0.02, which is 
different from that of the corresponding 3 + 1-dimensional case. In addition, the critical solution is also different, 
and, instead of having DSS, now has HSS. The above results were confirmed by independent numerical studies [10]. 
However, the exponent obtained by Husain and Olivier (HO), 7 ~ 0.81, is quite different from the one obtained by 
PC. It is not clear whether the difference is due to numerical errors or to some unknown physics. 

After the above numerical work, analytical studies of the same problem soon followed up [11-14]. In particular, 
Garfinkle found a class, say, S[n], of exact solutions to the Einstein-massless-scalar field equations and showed that 
in the strong field regime the n = 4 solution fits very well with the numerical critical solution found by PC. Lately, 
Garfinkle and Gundlach (GG) studied their linear perturbations and found that only the solution with n = 2 has one 
unstable mode, while the one with n = 4 has three [13]. According to Eq.(2), the corresponding exponent is given by 
7 = l/|fc|=4/3. Independently, Hirschmann, Wu and one of the present authors (HWW) systematically studied the 
problem, and found that the n = 4 solution indeed has only one unstable mode [14]. This difference actually comes 
from the use of different boundary conditions. As a matter of fact, in addition to the ones imposed by GG [13], HWW 
further required that no matter field should come out of the already formed black holes. This additional condition 
seems physically quite reasonable and has been widely used in the studies of black hole perturbations [15]. However, 
now the corresponding exponent is given by 7 = l/|fc| = 4, which is significantly different from the numerical ones. 
So far, no explanations about these differences have been worked out, yet. 

Self-similarity is usually divided into two classes, one is the discrete self-similarity mentioned above, and the other 
is the so-called kinematic self-similarity (KSS) [16], and sometimes it is also called continuous self-similarity (CSS). 
KSS or CSS is further classified into three different kinds, the zeroth, first and second. The kinematic self-similarity 
of the first kind is also called homothetic self-similarity, first introduced to General Relativity by Cahill and Taub in 
1971 [17]. In Statistical Mechanics, critical solutions with KSS of the second kind seem more generic than those of 
the first kind [2]. However, critical solutions with KSS of the second kind have not been found so far in gravitational 
collapse, and it would be very interesting to look for such solutions. We shall present in this work the study of the 
linear perturbations of the 2 + 1-dimensional circularly symmetric solution, obtained in a previous work [18], with 
kinematic self-similarity of the second kind and show that the background solution is not critical. 

In Section II we present the field equations with kinematic self-similarity of the second kind. In Section III we 
perturb linearly the field equations. In Section IV we present the solution of the linear perturbation equations. I 
Section V we apply the boundary conditions in the perturbed solutions and in Section VI we conclude our work. 



II. THE FIELD EQUATIONS WITH KINEMATIC SELF-SIMILARITY OF THE SECOND KIND 



The general metric for such spacetimes can be written in the form 

ds 2 = e 2*(r,t)^2 _ e 2*(r,t) rfr 2 _ ^2 

Then, the corresponding non-vanishing components of the Ricci tensor are 



„2(*-*) 



Rtt 

Rtr 



<J>, ( $ r - <Jr r + ^ + I ) + < TV( 



r * S S rS S 



= e 2(*-*) 



$r 



rb b r 



-2* 



!0 + * r -* r ) + i(* r -* r ) + ^' 



(3) 



(4) 
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where the indices t and r denote differentiation to the time coordinate and the radial coordinate, respectively. 
Then, introduce two self-similar variables r and x via the relations 



x = In 



r = -In(-t), 



or inversely, 

r = e ( M - T )/«, t = -e- T , 
where a is a dimensionless constant. For any given function / (t,r), we have 



(5) 



(6) 



/,t = -— 7 ( a .f.r + f.x) , f,r = -f.x, 

at r 

f.tr — 7 {ptf,rx ~t~ f.xx) i f.rr — o {f.xx ; 

f> tt = ^2 ("V.tt + 2af, TX + f iXX + a 2 f iT + af iX ) , 



(7) 



where the comma means differentiation. 

Substituting these equations into Eq.(4), we find that in terms of the self-similar variables, the Ricci tensor is given 

by 



2(S-*) r c 

Rtt = | - + -f ) + $ a 



f 



Y 2 t 2 



* TT + * T (1 + * T - $ T ) + ^ + ^(1 - 4> T 



+ a 
+ 



2* r:c + * T (** - $ x ) + * X (* T - $ T + 1) + 2^ + ^(1 - $ T ) - ^<P 2 

T ra <^ * xv x ^x) g g x 



Rtr 



1 



atr 

e 2(*-*) 



* r (l + |) + |(*.-D-^' 



Q 



S 



(8) 



-2$ 



f 



<I>, (*, - <!>,. -t I.) -<!>,,, r ( 1 + ~^(S XX + S X ) 



-2* r 



- (S x (1 + $ x - tf x ) + S xx ) + $ x - 



(9) 



where the indices r and x mean differentiation in respect to these variables. 



III. THE LINEAR PERTURBATION OF THE FIELD EQUATIONS 

Once we have the general expressions of R^ v in terms of r and x, let us consider the following perturbations, 

$(t, x) = $o(z) + e$i(x)e kT , 
*(r, .t) = * (a;) + e*i(x)e' CT , 
5(t,.t) = 5 (a;) + e5i(a ; )e feT , 

#r,a;) =^ o (t) + e0i(a:)e* T , (10) 
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where e is a very small real constant, the quantities with subscripts "1" denote perturbations, those with "0" denote 
the background self-similar solutions. 



The background solution is given by 



*o(a0 =0 
*o(*) = - 
S (x) 



\ax, 
2 



2-a 
00 (*) =2gln(-t), 



and the apparent horizon is given by 



rxff(t)= (2-a)(-t) 



1/2 



2/(2-a) 



(11) 



(12) 



where t/?o, -00 and s are integration constants, q = ±^75 an d a < 2 [18]. 

It is understood that there may be many perturbation modes for different values (possibly complex) of the constant 
k. The general perturbation will be the sum of these individual modes. Those modes with Re(k) > grow asr^oo 
and are referred to as unstable modes, while the ones with Re(k) < decay and are referred to as stable modes. By 
definition, critical solutions will have one and only one unstable mode. 

Substituting Eq.(10) into Eq.(8) and then we will have 



R^v — R^Lvi^-, x, e). 

Now considering is function of e only, then wc expand it in terms of e, 

R^(t, x, e) = — { JjW (x) + eR$ {x)e kr + O (e 2 ) } 



(13) 



(14) 



where R^i{x) is the part of the Ricci tensor corresponding background, and R$(x) the perturbation part, which 
is function of the background, &o(x), *&o, Sq(x) and the linear perturbations, $1(2;), Si(x). In the paper of 
Hirschmann, Wang & Wu [14] it was calculated them for the self-similar solutions of the first kind but in double 
null coordinates [cf. Eq.(65) given there]. 

To first order in e, it can be shown that the non-vanishing components of the Ricci tensor are given by 

R£\x) = e 2(*o-*o-*+-) | $o / (2$1 / _ ^ $1 ^ Q / + $1 « + 2 ( $1 - * l} « - + $ "] 
-V + 2 ($1 - *i) $o'S* ' + $ '5i' + *i'So 



1 



So Si ^ , 



(15) 



=2r 



+ — {$0' (afctf 1 + *i') + #o' + $/) - 2* ' (afctf 1 + */) - a 2 k 2 ^ 1 - 2afctf i' - *i" 



Q 



+ -L [$„' ( a A5i + Si') + So' (afc$i + $1') - a 2 fc 2 5i - 2akSi' - Si" - a 2 /eft - aSi' 
So 

1 



So 



r£\x) 



(So'Si$o' - S "Si - aSo'Si) 



Si 



aSo 



= 2(* -*o+t) 



a z 



-i (So' - $0%' - *o'S ' + So" - SoV) - $0' (afcSi + S/) - * 'Si' + 
00 

-S ' (afc*i + *i' + $i') + aA;Si' + Si" - * 'Si - S (afctf 1 + *i') + afcSi + Si'] 



(16) 



2 (*! - $!) (* ' 2 + * " + a<V - *„'*</ + 5-o'So') + 2afttfi*o' + 2tfo'*i' + <* 2 fc 2 *i + 
+ 2afc*i' + *i" + a 2 /j*x + a*i' - afc*i$ ' - *o'*i' - afc^o'^i - *o'*i' + 



4 



(akS *oSi + S *o'Si' + akS So'^i + S S '*i - S^q'Sq') 



+ 



s 2 (s-») 



-$i'*o' + *i" - $i' + 2$ '$i' - *o'*i' - tt (*o'5i ; + *i'5o' + Sb*i' - ft" - Si') + 



Si 



+ -i (* 'S ' - So" - So') 

^0 

= - e - 2 *° [(Si - 2S *i) ($ 'So' + $o'S - *o'S ' - *o'S + So" + So') + 

+S ($o'Si' + * 'Si + $i'S ' + *i'So - *o'Si' - * 'Si - - *i'5 + Si" + Si')] - (17) 

e 2(H_ir+,-* ) 

5 (Si - 2S $i) ($o'S ' - *o'S ' - So" - aSo') + 

or L 

+S (afc$ 'Si + $ 'Si' + a/c$iS ' + $i'S ' - akS^o - *o'Si' - afc*iS '+ 
— ^i'So' — a 2 k 2 S\ — 2akS\ — S\" — a 2 kS\ — aSi')] , 

where the prime denotes differentiation in respect to x. 
Once we have R$ (x) , we have to calculate the quantities 



Afiv = <t>,n<t>,v 



Substituting Eqs.(lO) into the above equations, we have 

V(w) = ^{A^(x)+eA^(x)e^+0(e 2 )}, 



(18) 



(19) 



where A^}(x) is the part of the background, and A$(x) the perturbation part, which is function of the background, 
(fro (t) , and the linear perturbation, (fry (x) , and given by 



2t 



A%\x) = -—\4q{ak<h+&) 



a 



a 



[2#i] 



4l\x) = - 

Ag>(x)=0 

where the dot means time differentiation. 

Once we have A$ (x) and R$ (x) , the linear perturbation equations are given by 



(20) 



(21) 



which in general are complicated. 

After we have the general linear perturbation equations (21), then we turn to consider the background solutions 
given by Eqs.(ll). By virtue of the simple form of the solutions and the fact (fro(x) = 0, Eqs.(21) can be solved in our 
case. 



a 2 fc$i + a$i' + a 2 k 2 ^x + 2ak^i + */' + 



„ (2 — a)e" x 



2k 2 + 2k + ^j a 2 Si + 2a (2k + 1) Si' + 2Si" 

— Aqa (akcfri + (fr[) ■ 



(22) 



*i" = 0, 



(23) 



(a - 2) [aft* i + + a*i' + - — [ a (2k + 1) Si + (2ak + a + 2) S/ + 2S/'] = -2#' l7 (24) 



2 k (2k - 1) *i + a (Ak - 1) + 2*i" + a (afc$i + *i') - a^- — [( 2a jfe + a ) Sj + 2S/] = 0, (25) 
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(2 - a) *!' + (2 - a) - 2$/' - ^-o^ 2 ' [aSi + (2 + a) + 25/'] - 0, (26) 
a (afc*i + - a (aft$i + $i') - (2-^) e2a: [a 2 /c (1 + 2ft) Si + a (1 + 4fc) 5/ + 25i"] - 0, (27) 
(2 - a) *i' - (2 - a) - ^-o^ 2 ' [aSi + (2 + a) Si' + 2S/'] = 0. (28) 



IV. THE SOLUTIONS OF THE LINEAR PERTURBATION EQUATIONS 

We will solve the system of the perturbed Eqs.(22)-(28). From Eq.(23) we have 

$i=aa; + b. (29) 

From Eqs.(26) and (28) we have 

(2 - a)*i = 0, (30) 
which solutions are a — 2 (which is out of range of our solution) or = 0. Thus, from Eq.(29) we have that 

$i=fe = constant. (31) 
Using Eq.(31) and summing the Eq.(26) and Eq.(28) we get 

P h ax 

*i = — + ( 2 + a ) S 'i + 2S "} ( 32 ) 
Using Eq.(31) and substituting Eq.(32) into Eq.(27) we get 



*i = b - -1 + 1 [a 2 fc(l + 2fc)Si + a(l + 4fc)5i + 



(33) 



Substituting Eq.(32) into Eq.(33) and differentiating it, we have 

ASi + BS[ + CS'{ + 4Si" + £e"5 aa: = 0, (34) 

where 

A= y(-8afc 3 + 4afc + a + 16fc 3 -4fc), (35) 

B = a(-8ak 2 +4ak + 3a + l6k 2 ), (36) 

C = 2(8a/c-a-4fc + 4), (37) 

£ = Aa 2 bk 2 . (38) 

Since our background solution with second kind self-similarity is identical to the solution with self-similarity of the 
first kind [18], we will study hereinafter only the case a = 1. Thus, 

^ = 4fc 3 + |, (39) 

B = 8k 2 + Ak + 3, (40) 
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C = 8fc + 6, 



(41) 



E = Abk 2 , (42) 

and the solution of equation (34) is given by 

Sl {x) = - b -ft- + Cie -^+ 2fc )* + C2e -i(fe+l+VA), + C3e -l 2 ( k+1 -VA )x _ (43) 
K 1 

where 

A = -fc(3fc-4). (44) 
In the next section we will apply the boundary conditions for two special cases: A > and A < 0. 

V. THE BOUNDARY CONDITIONS FOR THE PERTURBED SOLUTIONS 

We will apply the boundary conditions only at two regions of the spacetime: at the centre of the spacetime r = 
and at the event horizon tah given by equation (12), that furnishes 

tah = -t. (45) 

Thus, the metric at the apparent horizon is given by 

ds 2 AH = dt 2 -dr 2 -4(-t) 2 d9 2 . (46) 

It can be easily seen from this metric that the apparent horizon is singular, in this case, only at t — 0. Then the final 
state of the collapse is a marginally naked singularity. 

We would like to note that for the perturbed part of the metric (3) to represent circular symmetry, some physical 
and geometrical conditions needed to be imposed [19]. For gravitational collapse, we impose the following conditions 
at r = : 



(i) There must exist a symmetry axis, which can be expressed as 



X 



(47) 



as r — > , we have chosen the radial coordinate such that the axis is located at r = , and is the Killing 
vector with a close orbit, and given by ^" e jd a — dg. 

(ii) The spacetime near the symmetry axis is locally flat, which can be written as [20] 



AX 



1, (48) 



as r — > . Note that solutions failing to satisfy this condition sometimes are also acceptable. For example, when 
the left-hand side of the above equation approaches a finite constant, the singularity at r = may be related to 
a point-like particle [21]. 

(Hi) No closed timelikc curves (CTC's). In spacetimes with circular symmetry, CTC's can be easily introduced. To 
ensure their absence, we assume that the condition 



He) 

holds in the whole spacetime. 



< °> ( 49 ) 
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A. Case A > 



In this case we have from equation (33) that 

' [fc2 C3e ^(fc+V-M3fc-4)) + fc 2 C2e i*(fc-V-fc(3fc-4)) + ie xk kh _ 2c2fce ^(fc-V-fc(3fe-4)) 



1 



c 3 v/^fc(3fc - ^fcei^+V-fcC 3 *- 4 )) - 2c 3 fc e * a; (' £ +V- fe (3fe-4)) _ C2V /_fc( 3fc - 4 )fc e i*(W- fc ( 3fc - 4 )) 
c 3V /-fc(3fc - 4) e ^(fe+V- fe (3fe-4)) _ 2be xk + C2e ix(fe-V-fe(3fc-4)) + C2V /_ fc ( 3fc _ 4 )gix(fc-V-fc(3fc-4)) 



C3e 



U(fe+ A /-fc(3fc-4)) 



e + C4e 



(50) 



In order to apply the first boundary condition (47), we have to calculate the quantity \/X = rSi, which can be 
written as 



(-rt)i + c x r^ k (-t)^ +k + 



C2r -i( fe -i+VA) 



( _ t) i(Hi+#) +C3r -i(*-i-VA) ( _ t)i(t +i-v/i) i (51) 



i(fe-l-VA), 



i(Hl-VA) 



Since the limit of rSi must vanishes when r — > 0, all the exponents of r must be greater than zero. It is easily 
shown that some the exponents cannot satisfy this condition, then the first condition is not fulfilled. Thus, these 
perturbations are limited by the boundary conditions. 



B. Case A < 



In this case we have from equation (33) that 



*i(z) = - 



e 2 



xk 



fc-1 



-2c 2 fc 2 cos ( -<vA(3fc-4)z;) + 4c 2 fccos ( -^Jk{3k - 4)x ] - 3kbe^ xk 



+ 2c 2 ky / k(3k- 4) sin (]-^k{3k~A)x) - 2c 2 y / k(3k - A)sin ( \^/k{3k - 4)x ) + 26ei xfc 



2c 2 cos Qv/fc(3fc-4); 



+ c 4 e' 



(52) 



Then we have now two possibilities for the arbitrary constants c 2 and c 3 in order to get a real function S\{x): 
c 2 = c 3 and c 2 = -c 3 . 
For c 2 = c 3 we get 



Si(x) = - 



be , 2 X 
fc-1 



+ Cl e-5(l+2fe)- +2c2e -5(fe+l+)- 



COS 



1 



and rS\ is given by 



- +k V2c 2 r-^ k -^cos 



-A In 



H) 



(53) 



(54) 



Applying again the condition (47), we can see that all the exponents of r must be greater than zero only when 
k < 0, which admits only stable modes for the perturbation. For the second boundary condition (48), we have 
surveyed several sets of values of b, ci, c 2 , C4 and fc, and we have found at least one set that satisfies this condition: 
b = 0, ci = 0, c 2 = 1, C4 = 1, k = —1. In this case we get 



or 



lirrir^o rSi = -7 - 4V7, 
lim r ^ rSi = -4 + 4\/7. 



For c 2 = — c 3 we get 



Si(x) = - 



be z x 
k - 1 



Ci 



e -i(l+2fe)x _ 2 . C2e -i(fe+l+)x sin f i^/^A 



(55) 
(56) 

(57) 



Since this case is analogous to the case where c 2 = c 3 , we will not present it. 
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VI. CONCLUSIONS 



We have presented in this work the study of the linear perturbations of the 2 + 1-dimensional circularly symmetric 
solution, obtained by Chan, da Silva, Villas da Rocha & Wang [18], with kinematic self-similarity of the second 
kind. We have shown there that the solution is Kantowski-Sachs like [20] and it may be considered as representing a 
Friedmann-like cosmological model in a 2+1 dimensional spacetime. We have obtained in this paper an exact solution 
for the perturbed equations, which admits only stable modes, showing that our Friedmann-like background solution 
in 2+1 dimension is stable. This result is in agreement with the conclusions of the work of Hirschmann, Wang & Wu 
[14]. 
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